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Abstract 



We derive simple return models for several classes of bond portfolios. With 
only one or two risk factors our models are able to explain most of the return 
variations in portfolios of fixed rate government bonds, inflation linked govern- 
ment bonds and investment grade corporate bonds. The underlying risk factors 
have natural interpretations which make the models well suited for risk man- 
£h ' agement and portfolio design. 

1 Introduction 

Bond portfolios are subject to a large number of risk factors. Even in the case 
of default-free fixed rate bonds, the portfolio return depends in general on the 
' whole yield curve which is a high-dimensional object. Portfolio managers often 

, describe returns in terms of yield to maturity, which provides a one- dimensional 

^SJ ■ approximation of the expected return on a given portfolio. Such simple descrip- 

tions are easy to understand and to model, which makes them useful in risk 
management and asset allocation. 

This paper derives low-dimensional return formulas for several classes of 
bond portfolios, including fixed rate government bonds, index linked bonds and 
defaultable corporate bonds. The underlying risk factors in our models have 
natural interpretations which facilitates the assessment and modeling of the 
risks. For example, the returns on inflation linked bonds are quite accurately 
described in terms of the yield and the underlying consumer price index. Returns 
S-h ' on corporate bonds, in turn, are well approximated by the yield and the yield 

spread between corporate and government bonds. With only one- or two risk 
factors, our models are able to describe consistently over 97% of monthly return 
variations on government, inflation linked and corporate bonds over the past 
decades including the recent financial crisis. 

Our models are based on low-order Taylo r-approximations of t he log arithmic 
price of the portfolio. This is analogous to [Ch ance and Jordan! [l996| where a 
Taylor-approximation of the price with respect to yield and time was found to 
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give good approximations of return in the case of fixed-rate government bonds. 
Our model goes one step further by developing Taylor-approximations with re- 
spect to time, the yield as well as the outstanding coupon/principal payments. 
Variations in the outstanding payments are an essential return component e.g. 
in the case of index linked or corporate bonds. For inflation linked bonds, it 
amounts upto 10% increase in the explained historical return variations. Even 
in the case of default free fixed-rate bonds, our models improve on earlier ones. 
Compared to return formulas obtained by low-order approximations of the price 
itself, the logarithmic approach simplifies the formulas and improves the accu- 
racy of the approximations. 

2 Returns on bond portfolios 

Consider a bond portfolio whose coupon and principal payments are, as they 
are received, reinvested in bonds at current market prices. The portfolio's total 
return index Pt at time t then equals the market value of the fund's assets at all 
times. The yield to maturity, Yt of the portfolio is defined through the equation 

N 

P t = £e- y * (tn - i) Ct,n, (1) 

n=l 

where Ct,n denotes the fund's outstanding aggregate coupon and principal pay- 
ments payable at time t n > t. Whereas the default-free yield curve can be used 
to price an arbitrary portfolio of default free bonds, the yield to maturity is 
defined for a particular portfolio of (default-free or not) bonds. The yield to 
maturity of a portfolio is sometimes called the "internal rate of return" . The 
number JV of payment dates t n in flj can be allowed to be infinite but in that 
case our subsequent analysis requires that the outstanding payments Ct, n be 
bounded in n. An infinite N may be useful when modeling perpetual bonds. 

The definition of the outstanding payments is crucial in the definition of the 
yield. In the case of fixed rate bonds, there is little ambiguity: the outstanding 
payments are what the issuer has promised to deliver even though in the presence 
of default risk, the actual amount that will be received at time t n may be less 
than Ct,n- Default risk is reflected as a lower price and thus a higher yield. In 
the case of index linked bonds, the coupon payments are uncertain by definition 
and the outstanding payments in {TJ need to be forecast somehow. The simplest 
option is to "freeze" the underlying index to its current value when estimating 
Ct,n at time t. In the case of inflation linked bonds, the corresponding yield 
is known as the real yield. For now, it suffices to assume that the yield Yt is 
defined according to {TJ whatever the definition of outstanding payments may 
be. 

Equation |T]) expresses the bond price as a function of time t, the yield Y t 
and the sequence Ct = {Ct, n )t=i of payments outstanding at time t. Indeed, 
we have Pt = P(t,Yt,Ct), where P is the function defined for each t, Y and 
C = (GOSU by 

N 

P(t,Y,C) = J2 e ~ Y(tn ~ t)c - 

n=l 

Similarly, the log-return on the portfolio over a holding period [t, s] can be 
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expressed as 

AlnP « lnP(s,Y s ,C s ) — lnP(t,Y t ,C t ). 

The approximation is exact if there are no payments during the holding period, 
i.e. if ti > s. If s > t n , the nth payment has been collected and reinvested in 
the portfolio. An error may result if the reinvested payment appreciates during 
[tn,s] at a rate different from Y s . If the payments during [t,s] amount to a 
small fraction of the all outstanding payments (if length of the holding period 
is small compared to the maturity of the bonds in the portfolio), the resulting 
error could be expected to be small. 

Our reduced form models are based on Taylor-approximations of the log- 
return with respect to time, the yield and the outstanding payments. The first 
order approximation can be written as 

i ftp i ftp N i ap 

A ln p w T ~m {t ' Yt ' Ct) At + T W {t ' Yt ' Ct) AY + T da: {t ' Yt ' Ct)ACn 

n — l 

(2) 

where Ay = Ys — Yt and AC n = C a ,n — Ct,n- The first two terms are quite 
familiar. Indeed, we have 

1 BP 1 dP 

- — (t,Y t ,C t ) = Y t and — — (t,Y t ,C t ) = -D t , 

where 

t n=l 

is the Macaulay duration of the portfolio at time t. The first term of ([2]) rep- 
resents the pull to par effect, which means that bonds tend to appreciate when 
the maturity is approached. The second term gives the price sensitivity with 
respect to the yield. The effects of time and yie ld changes on returns off fi xed- 
rate government bonds have been studied e.g. in lChance and Jor dan 199(| who 
developed a Taylor series approximation for the price itself (as opposed to its 
logarithm) . 

The last term in © is the return component arising from changes in out- 
standing payments during the holding period [s,t]. Because of its multivariate 
character, the last term is difficult to model in general. However, when there 
are no portfolio updates during the holding period, there is often a single risk 
factor K a such that 

AC„ » K s C t ,„ (3) 

for all n = 1, . . . , N. For default-free fixed coupon bonds, ACt.n = so Q holds 
trivially with K s — 0. In the case of index linked bonds, on the other hand, 
K 3 equals the change in the underlying index over the holding period [t, s]; see 
Section [4] below. When © holds, the last term in ([2J reduces to 

N i ap i N 

n—l n—l 

and becomes 

AlnP ziYtAt-DtAY + Ks. (4) 
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The portfolio return over the holding period [t, s] is then approximated by 



« exp (Y t At - D t AY t + K.) - 1, 

where the right side is bounded from below by —1 just like actual portfolio 
return we are trying to approximate. This is why we linearized the logarithm of 
the price instead of the price itself as is often done. 

The logarithmic approach is supported also by the second order analysis. 
Denoting the convexity by 



C* = -^£(tn-t)V r ^-*)C t ,n, 



P 

n=l 

we get, by straightforward differentiation, that 

d 2 lnP , #lnP_ r n2 d 2 lnP _ e~^^ 

while all other second order derivatives of the logarithmic price are zero. If 
holds, the second order approximation of the log-return reduces (after simple 
algebraic manipulations) to 

AlnP « Y s At - D t AY + K s +hc t - D 2 )AY 2 (5) 

This differs from the first order formula Q only by the addition of the term 
quadratic in AY and in that, in the first term, Yt has been repl aced by Y s . This 
shoul d be compared with the more complicated return formula of lChance and Jordan! 
1996] who studied fixed-rate bonds through a second order approximation of the 
price itself (instead of its logarithm). The reduction of the second order terms 
in the quadratic approximation of the log return suggests that the logarithm of 
the price is well approximated already by the first order terms. 



3 Fixed rate government bonds 

In the case of fixed rate default-free bonds, we have ACt.n — so that © holds 
with K s ~ and ([SJ reduces to 

AlnPra Y s At- D t AY+hCt- D 2 t )AY 2 . (6) 

This is similar to the m odel studied e.g. in l llmanenl 1992] but there the time 
component was ignored. IChance and Jordan! Il99dl in corporated the time com- 
ponent but their model, like that of lllmanenl 1992l | . was based on a Taylor- 



approximation of the price instead of its logarithm. 



3.1 Empirical results 

We study the accuracy of the above models in explaining monthly returns on 
fixed rate government bonds. Our dataset covers end of month observations of 
the Barclays' market capitalisation weighted total return indices, durations and 
yields for France, Germany, Italy, United Kingdom, United States and the Euro 
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are4E The length of the time series for each country is given at the bottom of 
Table [T] 

We fit the following two models to the data 

Model 1: A In P a = c + Y s At - D t AY s + e a , 

Model 2: A In P a = c + Y s At - D t AY s + jAY t 2 + e a . 

The parameters c and 7 are estimated by ordinary least squares. The regression 
statistics provide us with diagnostic tools to evaluate the performance of the 
proposed models. By comparing the fits of the two models we can assess the 
significance of the quadratic term in explaining the returns. 

Table [1] displays the estimation results for the two model specifications. 
Model 1 already provides an almost perfect fit to the total return data with 
R 2 values ranging from 99.6% to 99.8%. The model fit has been consistent 
over time and the residuals have remained marginal even during the recent fi- 
nancial crisis. The R 2 values being close to 100% there is not much room for 
improvement when adding the second order term. The estimated coefficients 7 
in Model 2 a re generally significant at a 5% level but, consistently with previ- 
ous research [Chance and Jord an. 1996], the improvements in the R 2 statistics 
are marginal, less than 0.02% points in all the studied markets. The Partial- 
R 2 statistic is defined as the P 2 -statistic obtained by regressing the residual of 
Model 1 with the quadratic term. The quadratic term of Model 2 explains less 
than 5.11% of the residual variance of Model 1. The estimated constant c in 
Model 1 deviates from zero at 5% significance level, but with varying signs, This 
may indicate that the model has not entirely captured all the systematic return 
components. The addition of the quadratic term in Model 2 mitigates this effect 
to some extent. 



Table 1: Estimation results for fixed rate government bonds. 





FRA 


GER 


IT 


UK 


US 


EURO 


Model 1 














100 *c 


0.0147 


0.0083 


0.0084 


-0.0195 


-0.0202 


0.015 




(3.5842) 


(1.7215) 


(1.9339) 


(-4.8403) 


(-3.0658) 


(4.1467) 


R 2 


99.76% 


99.64% 


99.68% 


99.84% 


99.71% 


99.79% 


Model 2 














100 *c 


0.007 


0.0007 


-0.0001 


-0.0256 


-0.0277 


0.01 




(1.4255) 


(0.1186) 


(-0.0178) 


(-5.7509) 


(-3.6689) 


(2.2687) 


7 


24.4889 


24.1545 


32.5155 


4.8396 


11.1262 


17.031 




(2.7949) 


(2.2625) 


(2.5933) 


(3.0699) 


(1.9703) 


(1.9131) 


R 2 


99.77% 


99.65% 


99.70% 


99.85% 


99.72% 


99.79% 


Partial-P 2 


5.11% 


3.41% 


4.43% 


2.63% 


2.84% 


2.46% 


Data start 
Data end 


1997-12 
2010-03 


1997-12 
2010-03 


1997-12 
2010-03 


1980-12 
2010-03 


1998-12 
2010-03 


1997-12 
2010-03 



Notes. The table contains estimation results for the models: 

Model 1: AlnP, = c + Y s At - D t AY s + e s , 

Model 2: AlnP, = c+ Y s At - D t AY s + jAY t 2 + e a , 



Numbers in parentheses are the t-statistics for the estimated coefficients. 



Further information is available online at https://ecommerce.barcap.com/indices/index.dxml 
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Our results are consistent with lllmanenl [l992] who finds that the explanatory 
power of duration has increased over time. In his empirical study. Hlmanenl [l992|] 
found that the duration explained 80% to 90% of return variance of fixed rate 
government bonds during the 1980s. The nearly 100% i? 2 -values for Model 1 
during the past decade suggest that the trend has continued. It should be noted, 
however, that our results are not strictly comparable since our model explains 
the log-returns and it includes the time component. 

4 Inflation linked bonds 

The coupon and principal payments of most inflation linked bonds are tied to 
an underlying consumer price index so that the payment that will be received at 
time t n is Ct n , n ~ It„Co, n , where Co, n is the real payment at time t = and It„ 
is the value of the underlying consumer price index some time before t„. The 
specification of the indexation lag is specified in the contract and depends e.g. 
on the lag between time for which the index is computed and the time when 
its value is reported. We refer the reader to [Deacon et al.l . |2004| . Chapter 2] 
for a general overview of different cash flow structures on index linked bonds. 
The real yield of a portfolio of inflation linked bonds is defined by setting the 
outstanding payments according to the most recent value of the consumer price 
index It, i.e. 

Ct,n = ItCo,n 

in {T}. It follows that 

AC„ = IsC'o,n — ItCo,n 
= -K s AtCt,n, 

where ir s is the annualized rate of inflation over the period [t — S, s — 5] where S 
denotes the indexation lag. Assumption ([3]) thus holds with K 3 — ir 3 At, so the 
first-order approximation ([4]) can be written as 

AlnPw (Y s +iv s )At- D t AY. (7) 

4.1 Empirical results 

In our empirical study, we consider the two model specifications 

Model 1: AmP s = c + Y 3 At - D t AY 3 + e s , 
Model 2: AlnP s = c + (Y s + n s )At - D t AY s + e s , 

the first one of which ignores the third term in the Taylor-approximation ([2]). 
This allows us to evaluate the significance of the changes in the outstanding 
payments (the is'-term in ((4|) in explaining the portfolio returns. 

We use monthly observations of Barclays inflation linked government bond 
index data. The data consists of total return indices, yields and durations 
on portfolios of inflation linked government bonds for Canada, France, South 
Africa, Sweden, United Kingdom and United Statefl. The bonds' cash flows are 

2 Further information is available online at https://ecommerce.barcap.com/indices/index.dxml 
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linked to the evolution of country specific inflation indices (usually the general 
consumer price index) . The monthly time series of appropriate inflation indices 
were obtained from Eurostat and the national statistical authorities' websites. 
Since the indexation lag is not specified in the data, we select the lag 5 that gives 
the best fit to historical data. In all the studied countries, the data strongly 
supports a specific choice of the lag. The length of the time series for each 
country is given at the bottom of Table [2] 

Table [5] gives the regression statistics for the two model specifications in the 
six markets. The inclusion of the T^-term in Model 2 results in a substantial 
improvement over Model 1. Model 2 gives an extremely good fit to the total 
return data with R 2 values ranging from 97.6% to 99.8%. The lowest R 2 value 
for South Africa is mainly caused by one outlier observation caused by pricing 
distortions in the local conventional bond market in 2002 and the subsequent 
inversion of the breakeven inflation curve. 

The high values of Partial- J? 2 statistics mean that the K-term is able to ex- 
plain most of the residual variance of Model 1. Figure [1] depicts the unexplained 
residual returns for the two models. The reduction in the residuals between 
Model 2 and Model 1 is caused solely by the incorporation of the inflation term. 

The fact that the estimated constant terms in Model 2 are generally insignif- 
icant is another indication that the model is well specified. The country specific 
"optimal" lag lengths S are given on the last row of Table [2] 



Table 2: Regression statistics for inflation linked government bonds. 





CAN 


FRA 


SA 


SWE 


UK 


US 


Model 1 














100*c 


0.1642 


0.1328 


0.4935 


0.1114 


0.2073 


0.1966 


t-stat 


(6.201) 


(5.544) 


(9.367) 


(3.678) 


(7.611) 


(5.888) 


R 2 


97.48% 


95.46% 


78.6% 


91.29% 


96.06% 


94.41% 


Model 2 














100*c 


0.0046 


-0.0098 


-0.0069 


0.0087 


-0.0151 


-0.0004 


t-stat 


(0.628) 


(-0.839) 


(-0.389) 


(1.465) 


(-1.918) 


(-0.066) 


R 2 


99.80% 


98.91% 


97.58% 


99.66% 


99.67% 


99.82% 


Partial-i? 2 


92.22% 


76.02% 


88.69% 


96.14% 


91.43% 


96.78% 


Lag 


2 


2 


3 


2 


1 


3 


Data start 
Data end 


1996-12 
2010-03 


1998-09 
2010-03 


2000-03 
2010-03 


1996-12 
2010-03 


1996-01 
2010-03 


1997-03 
2010-03 



Notes. The table contains estimation results for the models: 

Model 1: A In P s = c + Y a At - D t AY s + e s , 

Model 2: AlnP s = c + (Y 3 + 7r s _ fc )Ai - D t AY a + e a , 



Numbers in parentheses are the t-statistics for the estimated coefficients. 



5 Corporate bonds 

It was found in lllmanen et all [19941 ] that duration's explanatory power in ex- 
plaining bond returns quickly decreases when moving from default-free bonds to 
corporate bonds. This section derives a simple two-factor model for the returns 
of well-diversified portfolios of non-callable corporate bonds. Our aim is to show 
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NcocnOT-wmTfm(D[^o)oi 

CTIOICDOOOOOOOOOO romo^cMmTmtD^coco 
CTiCJiCDOOOOOOOOOO mcnoooooooooo 



(a) Canada (b) France 




(c) South Africa (d) Sweden 




(e) United Kingdom (f) United States 

Figure 1: Residual returns (in percentages) of Model 1 (blue) and Model 2 (green) 
for inflation linked bonds. 

that the effect of defaults can be approximated by <(3j where K 3 can, in turn, 
be approximated by the yield spread between corporate and government bonds. 
We start by a mathematical justification based on familiar assumptions on the 
default process and on the diversification of the portfolio. The model will then 
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be validated by an empirical analysis of Merril-Lynch data on corporate bond 
portfolios. 

Because of defaults, the outstanding payments Ct.n of a portfolio of corporate 
bonds may change during the h olding period \t, s]. We will use the multiple 
defaults model described e.g. in ISchonbucherl |l998| ] and [Schonbucherl 120031 . 
Chapter 6]. The multiple defaults assumption can be justified by restructuring 
of defaulting issuers that often takes place in practice. Thus, we assume that 
each default event reduces the outstanding payments Ct, n to a fraction pCt, n 
of their pre-default values. The number p is called the recovery rate, which 
in the context of portfolio modeling represents the reduction in the portfolio's 
outstanding payments due to an individual default. A crucial assumption in 
justifying (|3}, is that each default event reduces the outstanding payments Ct, n 
for all n = 1, . . . , N by, approximately, the same recovery rate. 

The time of occurrence and the recovery rate associated with the ith default 
will be denoted by ti and pi, respectively. The sequence (ti,pi)°Z 1 then deter- 
mines the outstanding payments remaining after default events during a time 
interval [t, s] . Indeed, we have 

Is 

C s . n = J^J Pi Ct , n , 

i=l 

where I s = max{i j ti £ [t, s]} is the number of default events during [t,s]. 
Both the default times ti and the recovery rates pi are random in general. 
In order to allow dependencies between them, we will m odel (ti, Pi)^L i as a 
C ox process in [0, oo) x [0, 1] with intensity A t , p ; see e.g. ISchonbucherl |l998l ] 
or lGaspar and Slinkol |2008l |. That is, we assume that, conditionally on A, the 
number of default events with (ti, pi) in a measurable subset B of [0, oo) x [0, 1] 
is Poisson distributed with parameter 

Xt.pdtdp 

B 

and the number of events in disjoint subsets of [0, oo) x [0, 1] are independent . 
One could model defaults by a general point process (see e.g. [Kallenberd . 120021 
Chapter 12]), but for notational simplicity, we assume that the process can be 
expressed in terms of an intensity. Note that the default times follow a Cox 
process with i ntensity At = X t.pdp. The above m odel thus extends the one 
studied e.g. in lLandol [l998l ]; see ISchonbucherl |l998l | for further discussion and 
references. 

It follows (see the Appendix) that the amount of outstanding payments due 
at time s > t has conditional expectation 

E[C a , n I A] = Ct, n exp (- l u du^J , (8) 

where 

lu= (1 - p)\u, P dp 
Jo 

is the mean loss rate. Note that we can write U = At(l— pt), where At = J Q Xt.pdp 
and pt = Jq 1 pXt.pdp/Xt for At > and p t — other wise. This corresponds to th e 
defaultable term structure models proposed e.g. in lDuffie and Singleton! |l999l |. 
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Arguing like in Ijarrow et al,l |2005l ] (who studied defaultable zero coupon 
bonds with zero recovery), one can show that the default event risk can be 
diversified away in the sense that the outstanding payments C a ,n converge to 
E[C s ,n | A] when the number of issuers in the portfolio is increased; see the 
Appendix. This suggests that for a well-diversified portfolio, the effect of default 
events on outstanding payments can be approximated by 



C,„ « Ct, n exp I - J l u du ) (9) 

and thus that (|3]) holds with 

K a = exp ^— j l u du 

In order to estimate K 3 , we will assume that there exists a risk neutral measure Q 
under which the market prices of trade d securities are equal to the expectations 
of their discounted cash flows; see e.g. [Duffle and Singleton! . I2003I . Chapter 5]. 
The default-free short rate will be denoted by r. The amount of cash that will 
be received at time t n is equal to Ct„, n , the outstanding payments at time t n 
payable at time t n . Applying (J5J) with s = t n , we get 

Pt = E Q ex P (~ J t r « dM ) Ctn.n 

~ E Q exp (~ J rudu^jcxp^-J^ l u du^j Ct, n 
= E Q J2 exp ^- ^ " r c u du^j C\ 



y t,n , 



where r| = r t + h is the defaultable sh ort rate. This is an a logous to the reduced 
form pricing formulas obtained e.g. in ISchonbucherl [l998], Duffi e and Singleton! 
1999]. The instantaneous loss rate can thus be expressed as the short spread 
It = r c t -r t . 

Unfortunately, neither r£ nor rt is observable in practice. We will approxi- 
mate them by yields on short maturity corporate and government bond portfo- 
lios, respectively. This suggests the approximation 

K s ~ exp(-S s Ai) - 1, (10) 

where S 3 is the short spread at time s. When 5 S A is small, we can simplify this 
to 

K s w -S.At. (11) 

Figures [2] and |3] display historical monthly values of K s and S s At in percentages 
for European and US corporate bond markets during 1996/2-2010/1. The values 
of K s were obtained by solving K s from ((4]) for observed values of log-returns and 
yields on Merrill-Lynch investment grade corporate bond indices. In the figures, 
the spread S s is the difference between yields on portfolios of corporate and 
government bonds of 1-3 years to maturity. The values of S 3 At have remained 
below 0.7% which means that the relative error in using (JTTJ instead of (|10|l 
would have been less than 0.0025% points. 
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» 4 - - 0.6 



(a) Maturity 1-3 years 



(b) Maturity 3-5 years 





(c) Maturity 5-7 years 



(d) Maturity 7-10 years 



Figure 2: Historical evolution (in percentages) of S t At (blue) and residual returns of 
Model 1 representing default losses (green) during 1996/2-2010/1 for the US indices. 



Using in j4]) suggests the approximation 

AlnP^ (Y s - aS s )At - D t AY 



(12) 



for log-returns on corporate bond portfolios. We allow the constant a to deviate 
from one since even with the shortest maturity available (1-3 years in our nu- 
merical examples below), the spread S s is likely to overestimate the mean loss 
rate l 3 . Indeed, the yield on corporate bonds contains a premium for unexpected 
variations in the future loss rate. In addition to this systematic default risk there 
may be a premium for default event risk which remains when the diversification 
argument used to derive ([9)l is not exactly satisfied; see the Appendix. 



5.1 Empirical results 

We will again consider two different model specifications 
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(a) Maturity 1-3 years 



(b) Maturity 3-5 years 



1 



(c) Maturity 5-7 years 



(d) Maturity 7-10 years 



Figure 3: Historical evolution of St At (blue) and residual returns of Model 1 repre- 
senting default losses (green) during 1996/2-2010/1 for the European indices. 



Model 1: A In P a = c + Y s At - D t AY S + e 3 , 

Model 2: AlnP 3 = c + (Y a - aS s )At - D t AY 3 + e s , 

which differ by the term aS s , where S s is the yield spread on short maturity bond 
indices. The parameter a will be estimated from the data. We fit the models to 
monthly observations of Merrill-Lynch investment grade corporate bond indicex 
data. Our dataset covers the total return indices, yields and durations of market 
capitalization weighted investment grade corporate bond portfolios of maturities 
1-3, 3-5, 5-7 and 7-10 years from the US and European markets. As a proxy 
for the instantaneous loss rate we will use the yield spread of 1-3 year short 
maturity bonds. 

Table [3] displays the estimation results for the European and Table [4] for the 
US data under the two model specifications. Model 1, which ignores the effects 
of default losses on the outstanding payments, already explains the majority of 
the return variations with R 2 values ranging between 96.4% and 99.2%. The 
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Table 3: Return regressions for the EU corporate bond portfolios. 





Model 1 


Model 2 




1-3 Y 


3-5 Y 


5-7 Y 


7-10 Y 


1-3 Y 


3-5 Y 


5-7 Y 


7-10 Y 


100 * c 


-0.022 


-0.038 


-0.064 


-0.094 


0.013 


0.007 


0.028 


0.012 




(-3.431) 


(-5.014) 


(-4.357) 


(-5.754) 


(1.811) 


(0.818) 


(1.753) 


(0.675) 


a 










0.509 


0.633 


1.312 


1.502 












(7.424) 


(8.246) 


(8.921) 


(9.426) 


R 2 


96.38% 


98.63% 


97.5% 


98.21% 


97.33% 


99.05% 


98.35% 


98.87% 


Partial-P 2 










26.35% 


30.63% 


34.07% 


36.58% 


Data start 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


Data end 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 



Notes. The table contains estimation results for two different model specifications: 
Model 1: AlnP, = c+ Y s At - D t AY s + e s , 
Model 2: AlnP = (Y„ - aS B )At — D t AY s + e s , 



Numbers in parentheses are the i-statistics for the estimated coefficients. 



Table 4: Return regressions for the US corporate bond portfolios. 





Model 1 


Model 2 




1-3 Y 


3-5 Y 


5-7 Y 


7-10 Y 


1-3 Y 


3-5 Y 


5-7 Y 


7-10 Y 


100 * c 


-0.033 


-0.061 


-0.09 


-0.116 


0.02 


0.011 


-0.001 


-0.001 




(-3.506) 


(-5.265) 


(-5.424) 


(-8.145) 


(1.754) 


(0.747) 


(-0.035) 


(-0.098) 


a 










0.423 


0.571 


0.712 


0.915 












(6.462) 


(7.169) 


(6.005) 


(10.643) 


R 2 


97.81% 


98.51% 


98.35% 


99.19% 


98.27% 


98.88% 


98.66% 


99.53% 


Partial-i? 2 










21.33% 


25.02% 


18.97% 


42.38% 


Data start 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


1997-1 


Data end 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 


2010-1 



Notes. The table contains estimation results for two different model specifications: 
Model 1: AlnP s = c+ Y s At - D t AY s + e s , 
Model 2: A In P = (Y s - aS B ) At - D t AY S + e s , 



Numbers in parentheses are the t-statistics for the estimated coefficients. 



addition of the the spread term in Model 2 inmproves the model fit in all ma- 
turities and markets. The spread term is able to explain roughly 20-40% of the 
residual return variation of Model 1. In addition to improving the models' fit 
the addition of the spread terms also dilutes the significance of the estimated 
constant terms, which may be interpreted as a sign of an improved model spec- 
ification. These quantitative results are well supported by Figures and that 
display the evolution of unexplained residual returns of Model 1 together with 
the short maturity yield spread in the US and European markets, respectively. 

The estimates of a as well as the Partial-P 2 values tend to increase as a 
function of the index maturity. This can be explained, to a large extent, by 
the differences in the rating composition and seniority of the bonds included in 
the indices. Figure 3] displays the evolution of the rating distributions of the 
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US corporate bond indices by maturity bucket during the worst phases of the 
financial crises, when the default losses reached exceptionally high levels. The 
rating quality deteriorates with maturity which explains the increased default 
losses at the longer maturities. The composition of the US indices by instrument 
type are displayed in Figure[S] In line with the deterioration of the rating quality, 
the seniority of the bonds contained in the indices also declines as the maturity 
of the indices increases. The largest default losses were observed in the 5-7 year 
index which had the largest share of capital related debt instruments (under 
the Basel II regulation) and securitized exposures, both of which exhibited large 
losses during the financial crisis. On the other hand, the shorter end of the 
maturity spectrum contained the highest share of senior debt and thus exhibited 
the lowest default losses. Similar conclusions regarding the deterioration of the 
rating quality of the indices as a function of the maturity hold for the European 
indices as well but we omit the details. 



Illlllll 
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Figure 4: The composition of US corporate bond indices by rating classes AAA (dark 
blue), AA (light blue), A (yellow) and BBB (red) during 6/2009-5/2010. 
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(c) Maturity 5-7 years (d) Maturity 7-10 years 

Figure 5: The distribution of US corporate bond indices among senior debt (dark 
blue), subordinated debt (light blue), capital instruments (yellow) and securitized 
debt (red) during 6/2009-5/2010. 



The three terms in Model 2 explain majority of the return variation in the 
considered eight portfolios but, especially at the longer maturities, there remains 
residual spikes in the aftermath of the Lehman Brothers' collapse in September 
2008. Some idiosyncratic default event risk thus seems to remain, even though 
the indices contain bonds from hundreds of issuers. The systematic default risk 
on the other hand seems to be captured to a large extent by the spread term. 
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Appendix 



We pro ve jS} by us i ng th e notion of Laplace functional of a random measure; 
see e.g. [Kallenb crg. 2002, Chapter 12]. To this end, we write 



C s ,n — _i _i PiCt.n 
i = l 

= exp In p^j C t ,n 



= exp 



where f(u,p) = — X[t,s](u) hip and fi = H ere X[t,s] denotes the charac- 

teristic function of the interval [s,t] and 5t itPi is the Dirac measure that assigns 
unit mass to the point (ti,pi). Conditi onally on A , the p rocess (U, pOSi is 
Poisson (since it is a Cox process), so by [Kallenberd . 120021 Lemma 12.2(i)], 



E[C s , n I A] = exp 



J (1 - e- f{u ' p) )\ u , p dudp\ C t ,n 



exp y- J J (1 - p)\ u , p dudp ) C t) 
exp I - / l u du ) Ct, n , 



where 



I b 



I (l-p)X u ,pG 

Jo 



This gives (JSj . 

The approximation ([9]) can be justified by assuming as in ljarrow et al.l[2005l ]. 
that the portfolio is diversified among a large number M of issuers so that 

M 

c t ,n = y^tz/Ct^, 

where ui 1 is the relative weight and C\ n is the outstanding payments of the ith 
issuer, respectively. The defaults of each issuer are assumed to follow a Cox 
process with intensity A as described in Section [5] If, conditionally on A, the 
defaults of different issuers are independent, we get 

E{C s , n - E[C s>n I A]) 2 = EE[{C s ,n - E[C,, n \ A]) 2 | A] 

M 

= Ej2(u) 2 El(C%-E[Cl, n \\]) 2 \\] 

i=l 
M 

<E^(uCl n ) 2 , 

i=l 

where the inequality holds since since both C.'„. and E\C l_ n \ A] are bounded 
between zero and CI n . Assuming as in Ijarrow et all [2005],' that C\ n are uni- 
formly bounded in i and that the portfolio is well-diversified in the sense that 
J2iLi( wl ) 2 ~ * as M grows, we get that C a , n converges in mean square to 
E[C s , n \X}. 
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